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•^H . 1 INTRODUCTION 



ABSTRACT 

Chandra gratings spectra of a sample of 15 massive OB stars were analyzed under 
the basic assumption that the X-ray emission is produced in an ensemble of shocks 
formed in the winds driven by these objects. Shocks develop either as a result of 
radiation-driven instabilities or due to confinement of the wind by relatively strong 
magnetic field, and since they are radiative, a simple model of their X-ray emission was 
developed that allows a direct comparison with observations. According to our model, 
the shock structures (clumps, complete or fractional shells) eventually become 'cold' 
clouds in the X-ray sky of the star. As a result, it is expected that for large covering 
factors of the hot clumps, there is a high probability for X-ray absorption by the 'cold' 
clouds, resulting in blue-shifted spectral lines. Our analysis has revealed that such a 
correlation indeed exists for the considered sample of OB stars. As to the temperature 
characteristics of the X-ray emission plasma, the studied OB stars fall in two groups: 
(i) one with plasma temperature limited to ~0. 1-0.4 keV; (ii) the other wtih X-rays 
produced in plasmas at considerably higher temperatures. We argue that the two 
groups correspond to different mechanisms for the origin of X-rays: in radiative-driven 
instability shocks and in magnetically-confined wind shocks, respectively. 

Key words: stars: early-type - stars: winds - X-rays: stars - shock waves 



The present concept on the origin of X-rays in mas- 
sive OB stars posi ts that they are emitted by ho t gas 
heated by shocks (|Lucv fc White! 1 198(1 iLucvl Il982h . OB 
stars possess mas sive and fast winds driven by radia- 
tion pressure (e.g. ICastor. Abbott fc Kleinlll975l ) and sub- 
ject to instabilities (ra diation-driven instabilities, RDI, 
lOwocki fc Rvbickil Il984h which may give rise to the for- 
mation of strong shocks llOwocki. Castor fc Rvbickil 1 19881 : 



iFeldmeier. Puis fc Pauldrachl 1997 ). The resulting emission 
is expected to be relatively soft (plasma temperatures kT < 
1 keV), while the X-ray lines formed in the optically thick 
outflowing gas are predicted to be blue-shifted and asym- 
metric due to the different absorption o f the red and blue 
part of the line profiles l|lgnacel l200ll ; lOwocki fc Coher] 
l200lf) . 

Recently, it has been suggested that the presence of 
magnetic fields may be an important ingredient for the 
X-ray production mechan i sm in very young massive stars 
jBabel fc Montmerlei[l997l : lud-Doula fc Owockil 120021 ) . The 
main physical effect of an ordered magnetic field is to chan- 
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nel the stellar wind towards the equatorial plane where the 
flows from the opposite stellar hemispheres collide, leading 
to the formation of strong shocks. The basic feature of the 
so-called magnetically-confined wind shock model (MCWS) 
is that the X-ray emission should be much harder (plasma 
temperatures kT > 1 keV) than that originating from the 
RDI shocks, and that the X-ray line profiles should be nar- 
rower (line widths as small as j of the stellar wind velocity). 

With the launch of the Chandra and XMM-Newton ob- 
servatories, the quality of the X-ray data on OB stars has 
improved considerably, making it possible to test models of 
the origin of X-ray emission. Amongst the more than a dozen 
OB stars with gratings spectra of acceptable quality, only 
one object, £ Puppis, shows unambigu ous blueshifted spec- 
tral lines with a symmetric profiles l)Cassinelli et al.l 1200 ll : 
iKahn et alll200ll h On the other han d, there are two ob jects, 
9 1 Ori C (Schulz et al. 20 00, 2003: iGagne et alJl2005l ) and 
r Sco l|Cohen et all [20031 ). in which a clear manifestation 
of the MCWS effects is observed. While the scarcity of the 
MCWS-type objects might be accounted for by their youth 
and the expected decay of the magnetic field strength with 
the age of the massive star (|Schulz et alj|2003h . the standard 
RDI shock model needs further refinements to bring theory 
in accord with observations. 

The fact that most OB stars show no appreciable line 
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shifts and asymmetric profiles is an obvious indication that 
their winds are much more transparent in X-rays than origi- 
nally believed. This result, along with the small filling factors 
of the X-ray emitting plasma in OB stars as found in previ- 
ous (ROSAT) studies (|Kudritzki et alJll996h . lends support 
to a physical picture in which the stellar winds are clumpy 
and/or porous. Alternatively, if the winds are smooth and 
homogeneous, the mass-loss rates should have values appre- 
ciably s maller (a factor ^5 or more) than those presently ac - 
cepted jKramer. Cohen fc Owockll2003l ; ICoIien et alj|2006j ). 

Wind clumping and porosity effects on the X-ray emis- 
sion from OB stars, and specifically on the shape of the line 
profiles, have bee n explored in recent analytical and nu mer- 
ical models fe.g.. lFeldmeier. Oskinova fc Hamannl 20031; Os- 
kinov a, Feldmeier & Hamann 2004, 2006; Owocki fc Cohenl 
2006). ft has been shown that under given conditions (such 
as reduced mass loss and assumed 'typical' distance be- 
tween clumps) their inclusion may lead to results much more 
consistent with the observations. Given the complexity of 
these models and the fact that they are not yet fully self- 
consistent, it is important to gather empirical information 
about the physical conditions in the regions responsible for 
the X-ray emission to put further constraints on numerical 
models. 

Various diagnostics have been used in this re- 
spect: analysis of helium-like triplet ratios, global 
fits with discrete-temperature models, constructing 
a distribution of emission measures as function of 
temperature, based on fits to indivi dual-line fluxes 
and to the total X-ra y spectra JCassinelli et all 



20011: 



Kah n et all 120011: iWaldron fc Cassinellil 
Mille r etall 120021: I Cohen et all 120031: ISchulz et al.l 



2000; 



20041: Woidowski 



Leutenegger et all 



Sanz-Forcada. Franciqsini fc Pallavicinil 
fc Sc hulz 2004, 2005; lGagne et al.ll2005l : 
2006). These studies indicate that X-rays are produced close 
to the stellar surface (likely, in the wind acceleration zone), 
and that the corresponding hot plasma has a temperature 
stratification. The latter point reinforces the idea that the 
X-ray emission originates in an ensemble of shocks. 

Guided by this background, we have developed a simple 
model which bears all the basic charateristics of the X-ray 
production in wind shocks. The model is described in § [2] 
the data sample is given in § [3] the results are presented in 
§0] and discussed in §0 The conclusions close the paper. 



2 MODEL 

As in the case of the RDI and MCWS models, our basic 
assumption is that the X-ray emission of hot massive stars 
originates in shocks. Given the relatively high densities in 
the wind, the energy losses by the shock-heated plasma are 
considerable: thus, shocks should be radiative. This conclu- 
sion follows from simple estimates which show that the typ- 
ical cooling time of a parcel of gas at the postshock temper- 
ature and density is smaller than the typical dynamic time 
of the flow. Namely, for the shock position at a distance r 
from the star, the ratio of the cooling time (t c — ^p s h/Qc) 
to the dynamic time of the flow (td = r/v) is: 



tc 

td 



1.58 x lO" 2 ^ 6 



r 



and the ratio of the thickness of the radiative shock (i.e., 
the cooling length of a parcel of hot gas at the postshock 
temperature: l c = \v a ht c ) and shock 'radius' is: 



- = 3.59 x 10~ 3 T 3 V 
r 



r 
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V 1000 



Mr, 



(2) 



where T s h is the postshock temperature given in keV, viooo is 
the terminal stellar wind velocity (in units of 1000 km s _1 ), 
and Me is the stellar mass loss (in units of 10 -6 M© yr _1 ). 
For a strong shock, the relation between the postshock tem- 
perature (in keV) and the shock velocity (in units of 1000 km 
s _1 ) is given by T s h = 1.956pv 2 hock , and /i is the mean 
atomic weight. The relative number density of hydrogen is 
assumed xb = 0.9; thus, the relative electron number den- 
sity is x e = 1.1 for a fully ionized plasma. The gasdynamical 
quantities and the cooling function are described in § 12.11 

For typical mass-loss rates (Mq w 0.1 — 1) and wind 
velocities (wiooo ~ 1.5 — 2.5), the cooling time and cool- 
ing length of a shock developed in a wind are appreciably 
smaller than the corresponding typical characteristics of the 
flow (see eqs. pQ, [2]). This is true for postshock tempera- 
tures below 1 keV (i.e., shock velocities smaller than 1000 km 
s _1 for solar abundances) and shock locations from a few 
to several tens of stellar radii. Therefore, the assumption 
of a steady-state, plane-parallel, radiative shock is a good 
approximation for our analysis. This conclusion finds sup- 
port in numerical simulations of both the RDI and MCWS 
models. In fact, as a result of efficient cooling, the shocks 
'collapse' in geometrica lly thin shells and disk-like struc- 
tures, respectively (e.g. Feldmeier. Puis fc Pauldrachlll997l : 

it 
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(1) 



The description of our shock model and the ensemble 
of shocks is given below. The models were then used in the 
recent version (11.3.2) of the software package for analysis of 
X-ray spectra, XSPEC dArnaudlll996T ). A global-fit approach 
was adopted in our analysis for the following reasons: (i) the 
fit can automatically take into account the quasi-continuum 
due to numerous weak lines; (ii) by fitting the shape of 
the underlying continuum, the model places additional con- 
straints on the plasma temperature; (iii) the model can con- 
strain the column density of the X-ray absorbing gas; and 
(iv) it can yield estimates of relative element abundances. 
Finally, the X-ray emission is assumed to originate from a 
hot optically-thin plasma in coUisional ionization equilib- 
rium (CIE), as is the case for various types of astrophysical 
objec ts, including hot massive stars (e.g., IPaerels fc Kahnl 
l2006t) . 



2.1 Shock Model 

We consider a steady-state, plane-parallel, radiative shock 
moving in a gas with adiabatic index 7 = 5/3. The basic 
physical quantities of the flow (density or nucleon number 
density, velocity, and pressure) have their standard post- 
shock values for a strong shock: p s h = 4po (i s h = 4no)> 
v s h = Vo/4, and p s h — 3/ipov 2 , where the subscript '0' de- 
notes the preshock values (the gas velocity is given in the 
rest frame of the shock front). The mass and momentum con- 
servation equations that describe the downstream flow are: 
pv = pshVsh = povo (nv = n s hVsh = n v ), pv 2 + p = po^o, 
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Figure 1. Temperature profile in a steady-state radiative shock. 
Dotted lines denote the exact solution, while the dashed ones 
give the approximate solution under constant-pressure assump- 
tion. The four curves are for shock velocities of 500, 650, 750 and 
900 km s — 1 , as correspondingly labeled: 1, 2, 3, and 4. 



while the energy equation reads: 
-pv ■ 



d_ 



7 



7" 



1 J 



pv 
~2 



where £ is the distance downstream from the shock front, 
and Q c is the cooling due to the optically thin, hot plasma 
in CIE. 

In the steady-state case when cooling is efficient, the 
temperature decreases and the density increases downstream 
from the shock front. The cooling layer eventually 'col- 
lapses'; that is, the gas velocity asymptotically reaches zero. 
It follows (see the momentum equation) that the gas pres- 
sure does not experience large variariations in the postshock 
plasma (3/4 ^ p/poVQ ^ 1). Thus, we adopt a constant- 
pressure approximation (CPA) to describe the physical pa- 
rameters of the postshock plasma. As we will see below, this 
allows us to easily make use of the radiative shock model 
in the analysis of the X-ray emission. Under the CPA, the 
energy equation reads: 



5 dT 

-(1 + X e )nvk— = —XeXHU A(T) 



(3) 



where the cooling term is given explicitly, T is the plasma 
temperature, k is the Bolzmann constant, x e = n e /n is the 
relative electron number density, xh = ns/u is the rela- 
tive number density of hydrogen, and A(T) is the cooling 
function at CIE (e.g. iRavmond. Cox fc Smithll 19761 ). To fa- 
cilitate the calculations, we followed a standard approach 
and made use of a piece-wise power-law approximation to 
the cooling curve: A(T) = AqT q . In the high-temperature 
intervals (10 5 sC T < 4 x 10 7 , and T > 4 x 10 7 ) the adopted 
values of the constants are Ap = 7 x 10~ 19 ,3 x 10~ 27 and 



a — —0.6,0.5 (e.g., Mvasnikov, Zhekov fc Belovlll99ct ). 

Figure [1] shows a comparison between the exact and 
the CPA solutions for the radiative shock and illustrates 
that our choice (CPA) is acceptable for the purpose of the 
present study. Moreover, as we will see below, the CPA al- 
lows an easy incorporation of the radiative shock model in 



XSPEC which, unfortunately, is not the case for the exact 
solution of the problem. 

To determine the X-ray spectrum from a steady-state, 
plane-parallel, radiative shock (with 'surface' A), we must 
integrate the hot plasma emissivity over a range of plasma 
temperatures in the temperature-stratified postshock zone 
(Fig. [1]): 



Sp(E) = A 



e{E,T)x e x H n 2 d£ 



(4) 



where e(E, T) is the specific emissivity at given photon en- 
ergy E and plasma temperature T (its integral over the 
entire energy /spectral range gives the value of the cooling 
curve at T), f is the distance from the shock front, and T s h 
is the postshock temperature in keV. Then, using eq. ©, 
we can write: 



S P (E) 



f T sh 

F ah / e(E,T)T 1 - a d\nT 
Jo 



(5) 



where for numerical convenience we integrate over InT in- 
stead of T, and T s h is the normalization parameter. In this 
form, the X-ray spectrum from a radiative shock is well 
suited for modelling with XSPEC : namely, we simply have 
a known power-law distribution of emission measure at var- 
ious temperatures, and we adopt an optically-thin plasma 
model (model vapec in XSPEC ) for the plasma spectrum 
at temperature T (given by the quantity e[E, T]). Thus, our 
XSPEC shock model have the following parameters: the post- 
shock temperature T s h, chemical abundances and the nor- 
malization factor F s h, which is related to the total emission 
measure in the shock (some technical details are found in 
Appendix [A")) . 

2.2 Global Model 

Our global model assumes that X-rays are produced in 
radiative shocks which are randomly (and uniformly) dis- 
tributed in the stellar wind, most likely in the accelera- 
tion zone. Therefore, the spectral lines in the total (in- 
tegrated) X-ray spectrum from shocks are broadened by 
the bulk (wind) gas velocity, but are not shifted with- 
out including absorption. Qualitative physical consider- 
ations, also demonstr ated by numerical modelling (e.g., 
lOwocki fc Cohenll2001r ). show that the X-ray absorption ef- 
fects in the stellar wind result in blueshifted spectral lines 
with asymmetric line profiles. 

To describe this physical situation we have constructed 
a new XSPEC model, whose technical aspects are given 
in Appendix [X] The main features of the model are: (1) 
shocks are characterized by their postshock temperature 
and the distribution of the total emission measures in the 
shocks is determined from the Chebyshev polynomial algo- 
rithm, as is used i n the standard XSPEC model dopvmkl 
l|Lemen et al.lll989ri . We note that a similar approach was 
successfully adopted in fitting the X-ray emission from an 
ensemble of adiabatic shocks l|Zhekov et al.l [20061 ). (2) The 
basis vectors for X-ray emission from the shock ensemble 
are those from the radiative shock model incorporated in 
XSPEC (see § 12. ip . (3) All shocks have the same elemental 
abundances and share the same spectral line shifts. Thus, 
the basic model parameters are: the distribution of emission 
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Table 1. OB stars with Chandra Grating Spectra 



No. 


HD 


Name 


Type a 


Distance 6 


Observation 


Counts c 


Data d 










(pc) 


IDs 


(MEG) 




1 


150136 




03 


1300 


2569 


8569 


8 


2 


66811 


C PuP 


04 I 


450 


640 


17435 


1, 11 


3 


37022 


9 1 Ori c 


04-6p 


450 


3, 4 


28719; 13706 


3, 6, 7, 11 


4 


149757 


eoph 


09V 


140 


2571, 4367 


5915 




5 


36486 


5 Ori A 


09.5 II 


501 


639 


6116 


4, 11 


6 


37742 


C Ori A 


09.7 lb 


501 


610 


9220 


10, 11 


7 


37128 


e Ori 


BOIab: 


412 


3753 


6811 




8 


149438 


T SCO 


B0.2 V 


132 


638 


16838 


2, 11 


9 


206267A 




06.5 V 


800 


1888, 1889 


1597 


5, 11 


10 


47839 


15 Mon 


07 V 


767 


5401, 6247 


1631 




11 


57061 


r CMa 


09 II 


1480 


2525, 2526 


1470 


5, 11 


12 


37043 


i Ori 


09 III 


440 


599, 2420 


4934 


11 


13 


37468 


cr Ori AB 


09.5 V 


353 


3738 


1893 


9 


14 


111123 


P Cru 


B0.5 III 


110 


2575 


1830 


11 


15 


116658 


Spica 


Bl III-IV 


80 


4509 


3041 





a Spectral type is from Wojdowski & Schulz (2005) except for HD 150136, f Oph, e Ori, 15 Mon, 
a Ori AB, and Spica whose spectral type is from SIMBAD. 

b Values of the distance are taken from Wojdowski & Schulz (2005) except for HD 150136 (Herbst 
& H avlen 1977), and C Oph, e Ori, 15 Mon, a Ori, an d Spica (from the Jim Kaler's 'STARS' web 
site: |http:/ /www. astro. uiuc .edu/ ~kaler / sow / sowlist . h tml} . 

c The total number of counts in the Chandra MEG background-subtracted spectra. Thanks to 
the X-ray brightness of Ori C, the two observations were analyzed separately. The data for 
Spica were obtained with the Chandra LETG. 

d References to the a rticl es where the X-ray spec tra of the stars ar e dis cussed in detail: (1 ) 
ICassinelli et alj <200ll); (2) ICohen et alj d2003ll; (3) iGagne et al.l 120051); (41lMiller et alj ll2002t) ; 
(51 ISchulj [|2003r i; (6l lSchulz et alj teOOdl; (7llSchulz et alj d200ai: (8 1 ISkinner et alj d2005h; (91 
Skinn er et al. (2007, in preparation); (lOl lWaldron fc Cassinellil fcOOOh ; (lll lWoidowski fc Schulzl 
j2005h . 



measure of radiative shocks; chemical abundances, and the 
global line shift for the spectrum. 

The profile asymmetry is handled through two 
(half)Gaussians each approximating the blue- and red- 
halves of the line profile. The Gaussian widths can vary with 
wavelength (in the samer manner as in the XSPEC model 
gsmooth: AE oc E 13 ). Although the goal of our approach 
is just to obtain best fits to the line profiles over the entire 
spectrum, it is worth noting that the wavelength dependence 
of the line width represents the physical notion that faster 
shocks produce higher-temperature plasma whose emission 
is dominated by lines at shorter wavelengths. 



3 X-RAY DATA 

For our study, we selected a sample of OB stars with avail- 
able gratings observations in the Chandra data archive. We 
extracted first-order spectra for each object according to 
the procedure described in the Science Threads for Grat- 
ing Spectroscopy in the CIAO 3.3 Q data analysis software. 
The ancillary response functions for all spectra were gener- 
ated using the Chandra calibration database CALDB v3.2. 
Table Q] gives some basic information about the stars of the 
sample and the X-ray observations. Our analysis is based on 

1 Chandra Interactive Analysis of Observations (CIAO), 
htt p : / / cxc . harvard . edu/ciao/ 



the MEG spectra to take advantage of their higher sensitiv- 
ity compared to the HEG spectra. One object (Spica) was 
observed with the Chandra LETG. 

Depending on the quality of the data (the total num- 
ber of photons), the MEG spectra were rebinned to have 
between 15 and 30 counts per bin. The spectra near the 
intercombination and the forbidden lines in the helium-like 
triplets (Si XIII, Mg XI, Ne IX and O VII) were rebinned 
so that these two lines fall into one large bin. This technical 
approach is aimed at improving the quality of the fit when 
optically thin plasma models (as those available in XSPEC ) 
are used that do not take into account the effects of a strong 
external UV field. In the case of hot massive stars, such a 
UV field is present and it may alter the ratio of the two lines 
in the helium-like triplets, although it does not change their 
sum (the total emissivity/intensity). Thus, this technical ap- 
proach allows to preserve the piece of information available 
in the helium-like triplets that can contribute in constrain- 
ing some plasma characteristics (e.g., plasma temperature, 
abundances) which come out from the spectral fits. 



4 RESULTS 

The X-ray spectra (MEG) of the stars from our sample (see 
Table [1} were fitted with the global model that assumes a 
distribution of radiative shocks in the stellar wind (S I2.2[) . To 
determine the elemental abundances, we varied only those 
having strong emission lines in the observed spectum of a 
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given object. These usually include N, O, Ne, Mg, Si and 
Fe, while for some objects C and S were also varied. For 
the abun dance of all the other elem ents, the solar value was 
adopted (|Anders fc Grevessd Il989h . Finally, to derive the 
X-ray absor ption column density toward e ach object, we 
made use of iMorrison fc McCammonl (|l983T j cross sections 
(XSPEC model wabs). 

Table 2 presents some basic results from the fits to the 
spectra of the studied stars: the total \ 2 an< A the degrees of 
freedom for the fit (2nd column); the X-ray absorption col- 
umn density (3rd column); the individual elemental abun- 
dance relative to the solar value (columns 4-10); the observed 
and unabsorbed X-ray flux (last columns). And the corre- 
sponding distributions of emission measure of the radiative 
shocks in the stellar winds are shown in Figure [5] The qual- 
ity of the global fits in the vicinity of some strong emission 
lines is illustrated in Figure [3] A few properties are worth 
noting. 

There is an indication that the metal abundances are 
subsolar in the sample of OB stars with gratings spec- 
tra and this trend is well demonstrated by the iron abun- 
dance being a factor of 2-3 below the solar value. Nitro- 
gen might be considered an exception but we note that 
the poor quality of the data in the soft part of the X- 
ray spectra, where the nitrogen lines are found, prevents a 
firm conclusion. The values of the column density of the X- 
ray absorbrption material are consistent with those derived 
from the visual extinction to each obj ect, usin g data from 
iBerghofer. Schmitt fc Cassinellil |l996) and the Gorenstem] 
1 19751 ) conversion formula. The only object showing a clear 
sign of excess absorption with respect to its interstellar value 
is £ Pup which may suggest that wind absorption effects play 
an important role in this object. 

We see that most stars display the same distribution of 
emission measure of the radiative shocks, peaked at ~0.1- 
0.4 keV and with a very low level of high temperature tail 
extending above 1.0 keV. Thus, in most cases X-rays are pro- 
duced in low velocity shocks. This finding is consistent with 
the radiation-driven instabilities scenario. The two most no- 
table exceptions are 8 1 Ori C and r Sco. In the former case, 
the main peak is at energy ~4-5 keV with a weaker compo- 
nent at ~0.8 keV. The case of r Sco is more peculiar since 
there appear to be three peaks: one at very low energies, the 
main one at ~0.7 keV, and some emission above ~3 keV. 
One could argue that the hot tail in the distribution of emis- 
sion measure in these objects is due to the MCWS while the 
cool component comes from the RDI shocks. For example, 
the shape of the emission measure distribution in 6 1 Ori 
C is stable between the two data sets (taken at two differ- 
ent times) and only the total amount of the hot gas changed 
mostly in the high-temperature component (Fig. [2]). We note 
that this findi ng is in accord with the previous an alysis of 
Chandra data (|Schulz et a l. 2003; Gagne et al.ll2005r ). Such a 
different behaviour of the hot and cool components might be 
an indication of their different origin, as mentioned above. 
However, it is necesary to obtain more data also on r Sco, 
in order to establish whether such a behaviour is typical 
only for one object or has much more general validity. It 
should be noted that there is a third object in our sample 
of OB stars, £ Oph, whose distribution of radiative shocks 
is relatively 'hot' and it peaks at ~ 0.5 — 0.6 keV. 

The presence of considerably hotter plasma (therefore, 



fast velocity shocks), as definitely seen in 8 1 Ori C, r Sco 
and £ Oph, in combination with the relatively narrow lines 
observed in their spectra could be assumed as a clear man- 
ifestation of the ma gnetically-confined w ind effects. It was 
recently proposed bv lSchulz etail (|2003h that the efficiency 
of the MCWS mechanism correlates with the age of the mas- 
sive star. The relative youth of the first two objects (0.3 and 
~ 1 Myr, see the ir Table 4) and the small kinematic age o f 
C Oph (~ 1 Myr. iHoogerwerf, de Bruiine fc de Zeeuwll200ll ) 
indeed support this suggestion. 



5 DISCUSSION 

The radiative shock model introduced in this study allows us 
to deduce some important information on the global physical 
state of the X-ray emitting region of massive stars. As one 
can see from eq. (Al), the normalization factor of the flux, 
F s h, depends on the shock temperature (T a h), stellar wind 
parameters (Me, «looo), and distance to the object in kpc 

(d kpc ): T 3h = 2.444 8 ^P-, as 8 gives the effective 

surface of the shock. Thus, from the results of the fits to the 
X-ray spectra, namely, using the F 3 h value and knowing the 
distance to each star, we can deduce some useful information 
on the structure of the stellar winds. To do so, we introduce 

a new quantity, 8 ^ J^mo ) ' that is related to a specific shock, 
or to a group of shocks with the same velocity. We will refer 
to this quantity as the specific stellar 'cloudiness'. 

In our picture, the winds of hot massive stars are clumpy 
and the clumps are likely generated by radiation-driven 
instabilties, particularly efficient in the inner acceleration 
zone. The instabilities develop into strong shocks which are 
the source of X-rays. Since the shocks are radiative, the gas 
density in the cooling zone is much larger (probably by a 
factor of a few tens and more) than the local density of 
the smooth wind. Away from the acceleration zone, these 
shock structures (clumps, complete or fractional shells) ex- 
pand and cool down even further, but the density contrast 
in their 'tail' zone will remain appreciably high. Eventually, 
the clumps become 'cold' clouds in the X-ray sky of the 
hot massive star, and in a steady-state picture there will be 
a correlation between their characteristics: e.g., the higher 
the number of the hot clumps, the more numerous the 'cold' 
clouds; the higher the clumps' covering factor, the higher 
the chance for X-ray absorption by the 'cold' clouds. In such 
a picture, the X-rays emitted from a particular shock will 
likely experience some absorption by the cold gas in its cool- 
ing zone, and there is a chance for more X-ray absorption 
to take place in the 'cold' tails of the other shocks and the 
'cold' clouds as well. Although the exact situation depends 
on many specific details as the gas density in the clouds and 
their covering factor, it is anyway suggestive that the X-rays 
emitted from the 'far side' of the stellar wind (as seen by an 
observer) have a higher chance to suffer absorption than 
those emitted from the 'near side'. Therefore, the emergent 
X-ray emission might be altered appreciably, resulting in a 
blue-shift of the spectral lines. 

The specific stellar 'cloudiness' defined above is propor- 
tional to the quantities that determine the X-ray absorption 
efficiency of the cooling zones in the clumps and the cold 
clouds: namely, the geometrical extension of a cloud (the 



Table 2. Shock Fitting Results (MEG 1-st order spectra) 



Name 


X 2 /dof 




N (o) 








N 





















Ne 




Mg 




Si 


S 






Fe 




r X 


HD 150136 


146/335 


4 


2 ± 


6 






1 









38 


± 





13 





47 


± 


0.09 





66 ± 


10 


0.85 ± 0.11 


0.97 ± 0.31 





35 


± 


0.09 


3.9 (24.4) 


CPup 


491/433 


3 


7 ± 


2 


1 


69 


± 





33 





16 


± 





23 





38 


± 


0.03 





54 ± 


04 


0.85 ± 0.06 


1.36 ± 0.32 





34 


± 


0.03 


12.0 (152.) 


9 1 Ori C c 


776/1078 


4 


5 ± 


3 






1 









79 


± 





22 





68 


± 


0.07 





83 ± 


06 


1.05 ± 0.08 


1.11 ± 0.11 





59 


± 


0.04 


26.0 (53.0) 


COph 


111/223 





± 


9 


1 


43 


± 





63 





21 


± 





05 





51 


± 


0.09 





56 ± 


1 


0.65 ± 0.11 


1 





32 


± 


0.05 


5.7 (5.7) 


S Ori A 


179/241 





± 


9 





35 


± 





17 





27 


± 





07 





45 


± 


0.09 





51 ± 


11 


0.53 ± 0.12 


0.99 ± 2.26 





31 


± 


0.06 


10.7 (10.7) 


C Ori A 


209/234 





1 ± 


3 





23 


± 





09 





23 


± 





04 





44 


± 


0.07 





58 ± 


09 


0.63 ± 0.11 


1.48 ± 1.64 





36 


± 


0.06 


16.4 (18.5) 


t Ori 


130/256 





1 ± 


4 





33 


± 





12 





21 


± 





04 





39 


± 


0.06 





47 ± 


08 


0.58 ± 0.12 


1 





30 


± 


0.05 


11.5 (13.0) 


r Sco 


503/420 





± 


4 


1 


18 


± 





30 





29 


± 





04 





59 


± 


0.06 





63 ± 


06 


0.65 ± 0.06 


0.78 ± 0.24 





44 


± 


0.04 


16.8 (16.8) 


HD 206267A 


16/85 


2 


9 ± 1 


6 






1 









20 


± 





16 





27 


± 


0.16 





49 ± 


26 


0.61 ± 0.32 


1 





22 


± 


0.16 


1.1 (8.4) 


15 Mon 


25/83 





2 ± 1 


4 






1 









31 


± 





15 





59 


± 


0.33 





64 ± 


56 


0.76 ± 1.21 


1 





56 


± 


0.36 


2.5 (3.0) 


r CMa 


21/122 





± 1 


8 


1 


16 


± 


1 


99 





60 


± 





60 





71 


± 


0.60 





68 ± 


60 


1 


1 





49 


± 


0.44 


1.4 (1.4) 


i Ori 


134/201 





± 1 








79 


± 





36 





39 


± 





11 





64 


± 


0.15 





71 ± 


18 


1 


1 





56 


± 


0.14 


10.8 (10.8) 


a Ori 


47/141 





± 1 


8 






1 









44 


± 





26 





87 


± 


0.50 





94 ± 


57 


1 


1 





83 


± 


0.49 


2.3 (2.3) 


Cru 


46/132 





± 1 


9 





29 


± 





23 





14 


± 





08 





30 


± 


0.16 





21 ± 


17 


1 


1 





33 


± 


0.19 


5.0 (5.0) 


Spica d 


97/152 





± 


1 





04 


± 





07 





09 


± 





03 





23 


± 


0.08 





20 ± 


22 


0.16 ± 0.06 


0.42 ± 0.36 





22 


± 


0.09 


1.5 (1.5) 



Fit results from the XSPEC model: wabs(gsm) (gsm states for the global shock model as described in § 2.2). Derived for the each object emission measure 
~ distributions of radiative shocks are given in Fig. 2. The erros in the Table are la. All abundances are expressed as ratios to their solar values (Anders & 

^ Grevesse 1989). Those without associated errors were kept fixed at their solar value. 

a The X-ray column density is in units of 10 21 cm -2 . 
k-H b Flux units are 10 -12 ergs cm -2 s _1 . The observed X-ray flux (0.3 - 8 keV) is followed in parentheses by the unabsorbed value. 

c The results are from a simultaneous fit to the both observations of this object. Such a model shares the same abundances for the both spectra and the 
cj other model parameters were varied individually for each data set. The values for the column density and the X-ray flux in the Table are for Obsld 3 while 

g those for Obsld 4 are N H = 4.6 ± 0.4 and Fx = 19.0 (40.4). Abundances not given in the Table but also varied are: Ar = 1.42 ± 0.29 and Ca = 1.84 ± 0.47. 

d The model fit is to the 1-st order LETG spectrum. Carbon abundance is C = 0.22 ± 0.10. 

CO 
CO 
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Figure 2. Emission measure distribution of radiative shocks in the OB stars of our sample. Horizontal axes are the immediate 
postshock temperature in keV; vertical axes are the emission measure in units of 10 5e cm -3 . In the case of 9 1 Ori C, the results for 
Obsld 3 and Obsld 4 are shown by dotted and dashed lines, respectively. 



quantity 5) and the cold-gas density. The latter is propor- 
tional to the density of the smooth stellar wind, that is to 
the ratio of the mass loss to the wind velocity. Thus, one can 
use the total stellar 'cloudiness' in a given object, given 



2 Total stellar 'cloudiness': The distribution of emission mea- 



sure of radiative shocks is derived from the spectral fit of a given 
star. This yields the values of r s ^ for each temperature (T^). 
Then, based on eq. IjAlj l these values are used to derive the cor- 
responding specific 'cloudiness' (8) for each radiative shock in the 
distribution. Finally, the individual 8- values are added to give the 
total stellar 'cloudiness' for the studied object. 
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Figure 3. Portion of the background-subtracted spectra near some strong emission lines. The overlaid best-fit model is also shown (solid 
line). Horizontal axes are the observed wavelength (A); vertical axes are the flux density (photons s" 1 A~ 1 ). 



by the sum of all the specific values, as a measure of the 
overall efficiency of the cold X-ray absorbers in the wind. In 
such a case, a correlation between the total 'cloudiness' and 
the line shift of the spectral lines should be expected, in the 
sense that the larger the total 'cloudiness', the higher the 
amount of blue-shift of the X-ray spectrum of the star. 

Figure [4] presents the results of our estimate of the to- 
tal stellar 'cloudiness' vs. the global line shifts in the X-ray 



spectra of the massive stars of our sample. These results in- 
dicate that the expected correlation does exist: winds from 
stars with high values of the covering factor display a larger 
amount of blue-shift in the X-ray line profiles. Note that 
for many objects the derived from the fit errors can be ap- 
preciable. Unfortunately, the complexity of the model fits 
and the non-uniform quality of the data limit the accuracy 
on the derivation of the plotted qauntities. An additional 
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Figure 4. Total stellar 'cloudiness' vs. global line shift for the OB stars of the sample. The error bars (dotted lines) are la errors from 
the fit (the errors of the total stellar 'cloudiness' follow from the la errors for the total emission measure as derived in the fit, and the 
uncertainties of the distance to the studied object are not taken into account in the error propagation). The objects with good photon 
statistics (total number of counts in the spectrum higher than 5000) are marked by triangles. Data points within circles are binaries or 
multiple systems (as identified in SIMBAD). Stars are identified by numbers as in Table [T] The data point for 8 1 Ori C (number 3) is 
the average of the two observations. 



source of uncertainty comes from the 'cloudiness' parame- 
ter that depends on the accuracy of the cooling curve (see 
eq.[5}. However, the key point is that all the objects were an- 
alyzed in a uniform way. Thus, future improvements in the 
atomic data (i.e. cooling curve) might only shift the scale, 
but they will not affect the observed trend that we believe 
is rather robust. More importantly, as can be seen in Fig. [4] 
appreciable red shifts of the spectral lines are not observed. 

The observed correlation between the total stellar 
'cloudiness' and the spectral line shifts reveals additional 
details about the X-ray emitting zone in hot massive stars. 
For example, from Fig. U and for typical mass-loss rates 
(Ma « 0.1 — 1 ) and wind velocities (fiooo ~ 1.5 — 2.5), in 
most cases the 'cloudiness' parameter is small (8 ^ 1). In 
other words, dense clumps/clouds do not completely cover 
the X-ray 'sky' of massive stars. This also means that the 
RDI instabilities (and the MCWS as well) do not affect a 
big part of their stellar winds. 

Comparing the results for single and multiple sys- 
tems, we note that the latter show higher values of 
the 'cloudiness' parameter: this could result from addi- 
tional sources of X-ray e mission, such as colliding wind 
shocks in binary stars ( Luo. McCrav fc MacLowl 1990; 



Stevens. Blondin fe Pollock! 1 19921 : iMvasnikov fc Zhekovl 
19931 ). We have to keep in mind that the 'cloudiness' 
parameter is also a measure of the efficiency of the X-ray 
production. Therefore, if an extra source of X-rays exists 
in some object, then it will result in an effectively larger 



number of shocks representing the additional emission. This 
is the case for a binary system, where a higher value of that 
parameter is expected in comparison to a single star with 
the same spectral line shift. 

Finally, the fact that no appreciable redshift of the X- 
ray line profiles is observed in any object indicates that the 
distribution of the absorbers in OB stars likely has a spher- 
ical (or axial) symmetry. Future observations with a high 
spectral resolution and a good photon statistics will give 
us an opportunity to further test these conclusions and to 
study the X-ray gas kinematics in greater detail. They will 
also show whether the 'cloudiness' parameter for a given ob- 
ject changes with time which might be an important charac- 
teristic of the mechanism that drives the wind instabilities 
responsible for the X-ray emission from massive stars. 



6 CONCLUSIONS 

Using the Chandra public archive, we have analyzed the X- 
ray spectra of a sample of 15 massive OB stars. The basic 
assumption in this study is that the X-ray emission from 
such objects originates in shocks which develop in the stellar 
wind either as a result of radiation-driven instabilities or due 
to confinment of the wind by relatively strong magnetic field. 
The main results and conclusions are as follows: 

1. Based on the fast shock cooling in the winds of hot 
massive stars, we have developed a simple model of a steady- 
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state, plane-parallel radiative shock which has been incorpo- 
rated in the software package for analysis of X-ray spectra, 
XSPEC . 

2. The model is then used in a global analysis of the X- 
ray emission from the sample of OB stars with high quality 
spectra obtained by Chandra. We assume that the X-ray 
emission originates from an ensemble of shocks that are most 
efficient in the acceleration zone of the stellar wind. Far from 
the formation region the shock structures (clumps, complete 
or fractional shells) cool down and eventually become 'cold' 
clouds in the X-ray sky of the star. In case of high coverage, 
these clouds can effectively absorb the X-rays emitted by 
the underlying shocks. In such a case, the spectral lines in 
the X-ray spectrum of the massive star will be blue-shifted. 

3. Using our model, we define a new quantity, the stel- 
lar 'cloudiness', which represents the X-ray absorption effi- 
ciency by the clouds in the stellar wind. Our analysis reveals 
that despite the large intrinsic uncertainties, there is a cor- 
relation between the stellar 'cloudiness' and the global line 
blue-shifts in the observed X-ray spectra of the sample stars. 

4. For each star the model gives the distribution of emis- 
sion measure (DEM) of the radiative shocks responsible for 
the observed X-rays. The sample of OB stars with good spec- 
tra is distinguished in two groups: (i) one with the DEM 
peaked at ~0. 1-0.4 keV and with a very low level of high 
energy tail above 1.0 keV; (ii) one with the maximum of 
the DEM that falls at considerably higher temperatures. In 
the former case, the radiative-driven instability shocks are 
the likely mechanism for the production of X-rays, while in 
the latter (9 1 Ori C, r Sco and £ Oph) the X-ray emission 
originates in magnetically-confined wind shocks. 

5. The derivation of metal abundance indicates a sub- 
solar metallicity for all the stars of the sample, with iron a 
factor of 2-3 below the solar value. 
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APPENDIX A: THE MODEL 

The physical background of the radiative shock model was 
given in §[2] along with the basic features of our global model 
that assumes a distribution of radiative shocks in the stellar 
wind of massive OB stars. Some technical details about these 
two models are described in this Appendix. 

Shock Model: As discussed in § 12.11 the X-ray spectrum 
of a parcel of gas downstream from the shock front is given 
by the emission from an optically-thin plasma in CIE. Thus, 
the total spectrum of a radiative shock is an integral over 
a range of plasma temperatures, T ^ T s h, in the cooling 
zone (see eq. [5]), where T s h is the immediate postshock 
temperature for a strong shock. As seen from eq. J5), the 
emission measure in the cooling zone is a power-law function 
of temperature, T, and the spectrum at each value of T is 
assumed to be that given by the vapec model in XSPEC . 
The technical realization is based on the XSPEC subroutine 
sumdem with the APEC switch ON. For convenience, the 
numerical integration is over common logarithms (dlgT). 
The X-ray emission from plasma cooler than ~ 300, 000 K 
is not taken into account. 

Thus, our XSPEC model of the X-ray emission from a 
radiative shock has the following parameters: the postshock 
temperature T s h, chemical abundances, and the normaliza- 
tion factor T s h. As for all optically-thin plasma models in 
XSPEC , the model normalization is related to the emission 
measure of the hot plasma and in this case V s h gives the 
total EM in a radiative shock (for temperatures T T„h)- 

Global Model: Our global model assumes that the X- 
rays from massive OB stars originate in an ensemble of ra- 
diative slwcks_[£j32J : _We use Chebyshev polinomials algo- 
rithm (jLemen et aJjfl989h to describe the smoothed distri- 
bution of the emission measure in this ensemble of shocks. 
Since the basic parameter for shocks is the immediate post- 
shock temperature, this distribution is computed relative 
to T s h- Thus, the total EM in each shock of the ensemble 
is derived from the spectral fits. Note that the situation is 
similar in the XS PEC model c6pvmkl (which is based on 
ILemen et aIH i989') where the basic vectors are those giving 
the X-ray emission (spectrum) from optically-thin plasma 
in CIE. Here, the difference is that at each temperature T B h 
we have the X-ray spectrum from a radiative shock, as de- 
scribed in § 12. H and above. In the technical realization of the 
global shock model the total X-ray spectrum is derived by 
integrating the spectra from radiative shocks over the range 
of postshock temperatures. For convenience, this integration 
is over common logarithms (d\gT s h). 

It should be noted that there is one important qualita- 
tive difference between the use of the global shock model or 
the model with distribution of emission measure of optically- 
thin plasma (e.g., c6pvmkl in XSPEC ) for interpreting the 
observed X-ray spectra. In the former case, the EM at a 
given temperature, T s h, in the distribution is the total EM 
(including plasma with T ^ T s h) in a radiative shock with 



velocity corresponding to this postshock temperature. In the 
latter case, the EM is that of an isothermal hot gas at that 
given temperature in the distribution. If an object is stud- 
ied whose X-ray emission originates in radiative shocks, as 
is the case of massive OB stars, and a model like c6pvmkl in 
XSPEC is adopted, it is not quite correct to assign a shock 
velocity to each temperature value in the hot-gas distribu- 
tion. This is because the hot plasma in a radiative shock 
is not isothermal (§ 12.111 . Thus, the amount of hot gas at 
each temperature in the c6pvmkl distribution is the total 
contribution from all radiative shocks having some plasma 
at that temperature. This limitation is not present in our 
global shock model discussed in this work. 

Additional advantage of the radiative shock model is 
that its normalization factor, T s h, which gives the total emis- 
sion measure in the shock, can be expressed through some of 
the basic shock parameters (as done in deriving eq. [5] using 
eq. [3]). In the simplified picture of spherically-symmetric 
stellar winds, it is justified to assume that the shocks move 
radially, and their surface can be presented as a fraction, 5, 
of the 'local' sphere: A = S4tyR 2 . Thus, the normalization 
parameter, F s h, of our XSPEC model can be written: 

5A(l + x e )n ah v ah k IP" 14 
sh ~ 2A 47rd 2 

= 2.444 S (J^)^ (Al) 

where dk pc is the distance to the object in kpc, and all other 
quantities are defined in § [2] The quantity (10~ 14 /47rd 2 ) 
comes from the XSPEC normalization factor of the flux for 
optically-thin plasma models. Note that T s h does not de- 
pend on the shock radius since in the case under consid- 
eration the gas density, n s h, goes as R~ 2 . Again, the rela- 
tive number density of hydrogen is assumed xh = 0.9 and 
x e = 1.1 for a fully ionized plasma. For simplicity, in deriving 
eq. (IA1|) we replaced the local wind velocity with the termi- 
nal speed which means that the effective surface of the shock 
(presented by ft) has the 'maximum' value. We note that the 
quantity 8 may have values larger than unity since it repre- 
sents the total area of the shocks with a particular velocity 
within the entire ensemble of shocks. Thus, once the value 
of the flux normalization parameter, T s h, is derived from a 
spectral fit, it can be further used to gain insight into the 
physical picture, provided some of the stellar parameters are 
already known from other studies. 



